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Abstract 

A symmetric version of the well-known SPM model for sandpiles is 
introduced. We prove that the new model has fixed point dynamics. 
Although there might be several fixed points, a precise description of the 
fixed points is given. Moreover, we provide a simple closed formula for 
counting the number of fixed points originated by initial conditions made 
of a single column of grains. 

Keywords: SOC systems; sandpiles; fixed point dynamics; discrete dynamical 
systems. 

1 Introduction 

Self-Organized Criticality (SOC) is a very common phenomenon which can 
be observed in Nature. It concerns, for example, sandpiles formation, snow 
avalanches and so on p]. 

Practically speaking, it can be described as follows. Consider an evolving 
system. After a while, the system reaches a critical state. Any further move 
from this critical state will cause a deep spontaneous reorganization of the whole 
system. No external parameter can be tuned to control this reorganization. 
Thereafter, the system starts evolving to another critical state and so on. 

Sandpiles are a very useful model to illustrate SOC systems. Indeed, consider 
toppling sand grains on a table, one by one. Little by little a sandpile will start 
growing and growing until the slope reaches a critical value. At this moment, 
any further addition of a single sand grain will cause cascades of grains and 
deep reorganization of the whole pile. Afterwards the sandpile restarts growing 
to another critical state and so on. 

A formal model for sandpiles, called SPM, has been introduced in [21I3IHI- 
The sandpile is represented by a sequence of "columns". Each column contains 
a certain number of sand grains. The evolution is based on a local interaction 
rule (see Section ISJ: a sand grain falls from a column A to its right neighbor B 
if A contains at least two grains more than B; otherwise there is no movement. 
The SPM model has been widely studied CIl El CHI El ■ In particular, it 
has been proved that it has fixed point dynamics and a closed formula has been 
given to calculate precisely the length of the transient to the fixed point [^. 
Moreover, a precise description of the fixed point has been given [2l- 
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All these results are very interesting but they have two main drawbacks. 
First, they lack generality; indeed, the fixed point results are always obtained 
starting from very special initial sandpiles (just one column). In jH we 
tried to solve this problem by giving a fast algorithm for finding the fixed point 
starting from any possible initial condition. Second, the model lacks symmetry; 
in fact, grains either stay or move to the right only. Remark that in Nature, 
sandpiles evolve absolutely in a symmetrical manner. 

In this paper we introduce SSPM: a symmetric version of SPM. The new 
model follows the rules of SPM but it appHes them in both directions. For 
technical reasons that will be clearer later, we allow only one grain to move per 
time step. 

We prove that SSPM has fixed point dynamics. This is not a great surprise. 
To validate the new model, one should give a precise description of these fixed 
points and compare their "shape" with those of sandpiles in Nature. 

To this extent we use a formal construct which allows a better description 
of the dynamics: orbit graphs. They are directed graphs of the relation "being 
son of". In Section 18.21 the precise structure of their vertices is given (under 
the condition of considering initial configurations made by a single column): 
a configuration belongs to some orbit graph if and only if it admits a crazed 
LR-decomposition (see Section IT2I1 . 

Practically speaking, a configuration admits a crazed LR-decomposition if it 
can be decomposed into an increasing part L and a decreasing part R and both 
in L and in R any two plateaus (i.e. consecutive columns of identical height) are 
separated by at least a "cliff" {i.e. consecutive columns with height difference 
strictly greater than 1). 

The special structure of the vertices allows a very useful description of the 
fixed points: they are configurations which admit a crazed LR-decomposition 
without cliffs. 

Finally, using this characterization of the "shape" of fixed points we provide 
a closed formula which computes the number of fixed points originated from 
the initial configuration (n) (a single column containing n grains) . The surprise 
is that the formula is [-y/nj- Unfortunately, we have no practical or "visual" 
explanation for such a formula. 

2 The SPM model 

A sandpile is a finite sequence of integers (ci, . . . , c^); fc G N is the length of the 
pile. Sometimes a sandpile is also called a configuration. Let C = UfceNC*^^)*^ 
be the set of all configurations. 

Given a sandpile (ci, . . . , Cfe), the integer n — Si=i '^i is the number of 
grains of the pile. Given a configuration (ci, . . . , Cfc), a subsequence . . . , Cj 
(with 1 < « < j < fc) is a plateau if Ch = Ch+i for i<h<j; s = i — j + 1 is the 
length of the plateau and p = Ci its height. A subsequence Ci,Ci+i is a cliff if 
Ci - Ci+i > 2. 

In the sequel, each sandpile (ci, . . . ,Ck) will be conveniently represented on 
a two dimensional grid where q is the grain content of column i. 

A sandpile system is a finite set of rules that tell how the sandpile is up- 
dated. SPM jH] {Sand Pile Model) is the most known and the most simple 
sandpile system. All initial configurations contain n grains in the first column 
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and nothing elsewhere i.e. they are of type (n). It consists in only one local rule 
which moves a grain to the right whenever there is a cliff (see Figure QJ. 




Figure 1: local rule of SPM. 

Formally, for any configuration c, if there exists i e N such that c — q+i > 2, 
then c evolves to c' according the following relations: 

f 4 = Q-l 
I = Q+i + 1 . 

This process is iterated until the rule cannot be applied anymore. We say 
that a fixed point is reached. 

Along the evolution of the pile, the rule may be applicable at different places 
in the configuration. To illustrate this, we represent the set of reachable config- 
urations (starting from a single column) on an oriented graph where the vertices 
are the configurations. There is an edge between two configurations and 
when c? can be obtained by applying the local rule somewhere in (see Figure[3 
for an example, starting from a single column with 8 grains). This is called the 
orhit graph of the initial configuration c, denoted by Qc- 
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Figure 2: Q(%), orbit graph of a single pile with 8 grains for SPM. 

The following theorem proves that the fixed point is unique, independently 
of the order of application of the local rule. 

Theorem 1 ([01) • For any integer n, for SPM is a lattice and is finite. 

The following lemma characterizes the elements of the lattice. 

Lemma 2 (0). Consider a configuration c and let n be its number of grains. 
Then, c e G(n} for SPM if and only if it is decreasing and between any two 
plateaus of c there is at least a cliff. 

Remark 1. Consider a configuration c, and assume that c contains a plateau 
of length 3. Such a plateau can be seen as two consecutive plateaus of length 2. 
Thus, by Lemma\^ c does not belong to any orbit graph. 
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From Lemma 12 it is easy to see that a fixed point 11 is a decreasing config- 
uration with no chffs and at most one plateau. Therefore for any n € N, we can 
describe the fixed point 11 of (n) by 



n = 



{P,P- 1, • • • , 1) if g = , 

{p,p - 1, . . . ,q + l,q,q,q ~ I, . . . ,1) Otherwise, 

where (p, q) is the unique decomposition of n in its integer sum: 

, p-{p+l) 
n = q + 2_^i^q+ . 



3 The symmetric model 

In this section we extend SPM to SSPM (Symmetric SPM) according to the 
following guidelines: 

• a grain can move either to the left or to the right, if the difference is more 
than 2; 

• when a grain can move only in one direction, it follows the SPM rule 
(right) or its symmetric (left). 

For all configurations c = (ci, . . . , c^-), the following local rules formalize the 
above requirements: 




K(ci,...,Cfe) 



otherwise, 



(ci, . . . ,Ci_i+l,Ci-l, . . . ,Cfc) if 17^1 , 
(1, ci — 1, . . . , Cfc) otherwise. 



Let i5[(c) denote the difference between the grain content of column i and 
the one of column i + 1 of c; define Sl{c) = Ck- Similarly, (5,'(c) denotes the 
difference between the grain content of column i and the one of column i — 1 of 
c with S[{c) = ci. 

Notation. For a, 6 e N with a < b, let [a, b] denote the set of integers between 
a and b. 

From the local rule we can define a next step rule / : C i-^ *P(C) as follows 

fic) = {v:{c)\S:{c)>2,^e[l,k]}u{V,'ic)\Sl{c)>2,^e[l,k]} . 

Finally, using the next step rule, one can define the global rule which de- 
scribes the evolution of the system from time step t to time step t + 1: 

yse^ic), /(5) = U/(c) . 

ces 

When no local rule is applicable to c, i.e. /({c}) = 0, we say that c is a fixed 
point of SSPM. For n S N , let /" denote the 7i-th composition of / with itself. 

The notion of orbit graph can be naturally extended to the symmetric case 
by using the functions and . In the sequel, when speaking of orbit graph, 
we will always mean the orbit graph w.r.t. the SSPM model. 
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3.1 Fixed point dynamics 

In this section we prove that SSPM has fixed points dynamics. This result 
is obtained by using a "potential energy function" and by showing that this 
function is positive and non-increasing. 

Given a configuration c = (ci, . . . , Cfc), the energy of a column Ci {i G [1, fc]) 
is defined as follows 

i=i 

Therefore, the total energy of a configuration c = (ci, . . . , Cfe) is naturally defined 
as 

k 

Lemma 3. Consider a configuration c ~ (ci,...,Cfc) with n grains. Then it 
holds that E(c) < E{{n)); equality holds if and only if c = (n). 

Proof Remark that n = I]i=i'^«- Define h{i) = Y.'jIZT+c2+-+c,-i+i^- Then, 
E{{n)) = Y^^=ij can be rewritten as E{{n)) — Y^\^ih{i). Note that h{i) > 
e(ci) for any i € [1, fc]; equality holds if and only if i = 1. □ 

The function E can be naturally extended to work on set of configurations 
as follows 

V5 e «|?(C), E{S) = max{£;(c), c G 5} , 

with E{^) = 0. 

The following lemma is straightforward from the definition of the energy 
function. 

Lemma 4. For any set of configurations S ^ %, E{f(S)) < E{S). 

The following simple proposition describes the general structure of the orbit 
graph. 

Proposition 5. For any initial configuration c, Gc is finite, contains at least a 
fixed point but no cycles. 

Proof. If /(c) =0 i.e. c is a fixed point, then we are done. Assume that c 
is not a fixed point. Remark that its energy is finite. By Lemma ^ it holds 
that E{{c}) > E{fi{c})) and E{f{{c})) > E{f*+^{{c})) for t > 1 (unless 
/*({c}) ~ 0)- Since E is a positive function, there must exist /i e N such that 
ff^ilc}) = 0. Then, /''"^({c}) contains a fixed point. If /''({c}) = 0, then the 
orbit graph Gc is finite, since |/*({c})| is finite for any t G N. 

Finally, there are no cycles in Gc otherwise the elements of the cycle would 
contradict Lemma 21 □ 

The following corollary is given only to further stress the result of Proposi- 
tion 

Corollary 6. SSPM has fixed point dynamics. 
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Figure 3: ^(5) , orbit graph of a single pile with 5 grains for SSPM. Remark that 
there are two distinct fixed points. 



Corollary El says that independently of the order of application of local rules 
both with respect to type of rule and to the application site, SSPM evolves 
towards a fixed point. The problem is that this fixed point might not be unique. 
Figure gives an example of this fact. 

Despite the non-uniqueness, in the next section we give a precise character- 
ization of the structure of the fixed points. This characterization is essentially 
deduced from the properties of the vertices of the orbit graphs. 



3.2 Orbit graphs 

In [3, the authors precisely described the structure of the orbit graph of SPM 
when started on initial condition (n). They proved that it is the graph of a 
lattice. As a consequence, they deduced the uniqueness of the fixed point for 
SPM. 

We have already seen that in the SSPM case, the dynamics is of fixed point 
type, but the fixed point might not be unique. Hence, it is clear that the orbit 
graph of SSPM is no more the graph of a lattice. In this section, we detail the 
overall structure of the vertices of these graphs. 

A configuration c ~ (ci, C2, . . . , Cfc) is LR- decomposable if it can be divided 
into two zones: L{c) ~ [1, t], i?(c) = [t + 1, fc] such that 

1. Vi G L{c),i 7^ t,Ci < Ci+i i.e. L{c) is non-decreasing; 

2. Vi e R{c),i ^k,Ci> q+i i.e. R{c) is non-increasing. 



Figure 4(a) give an example of LR-decomposition. For any configuration c, 

let T{c) = {i e [1, k],\/j G [1, A;], Ci > Cj}. In the sequel, T(c) is called the top of 



c, see Figure 4(b) 



Given a configuration c = (ci, C2, . . . , c^), a set of consecutive indexes / C 
[l,k] is crazed if any two plateaus in / are separated by at least a cliff. A 
configuration c has a crazed LR-decomposition if it admits a LR-decomposition 
in which both R{c) and L{c) are crazed. 

A configuration might have several different LR-decompositions. The fol- 
lowing propositions tell which of them we are interested in. The proof of Propo- 
sition El will be made progressively, using several technical lemmas. 

Lemma 7. Consider n G N and c G G{n} ■ Then c is LR- decomposable. 



L(c) 



R(c) 



T(c) 



Q 




LR- 



(b) The top of c. 



Figure 4: decomposition of a configuration c. 



Proof. The thesis is trivially true for the initial configuration (n). Now assume 
that c € S(n). Let L{c) = [l,t],i?(c) = [i + l,k] be a LR-decomposition of c. 
Consider d £ f {c) and assume that d = ^^''(c). We have three cases: 

1. i = then L{d) i(c) \ {t} and R{d) = i?(c) U {t}. L{d) is non- 
decreasing since L{d) C L(c). is non-increasing since dt > dt+i and 
i?(c) C 

2. i e [t + - 1]; then = L(c) and = i?(c). Of course L{d) is 
non-decreasing and R{d) is non-increasing. 

3. I = fe; then L{d) = L{c) and R{d) = R{c) U {fc + 1}. Of course L{d) is 
non-decreasing. R{d) is non-increasing since R{c) C R[d) and dfe > 1. 

The proof is similar \{ d~Vl (c) . □ 

Lemma 8. Consider n e N and c G iei T(c) 6e i/ie iop o/ c. j4nj/ 

decomposition of c is such that both L{c)\T(c) and R{c)\T{c) have no plateaus 
of size strictly greater than 2. 

Proof. Consider ti G N and c € G(n) ■ Remark that if c = (n) then the thesis 
is true. Now, assume that c ^ (n), then c should have an ancestor in G(n)- 
We prove the thesis for R{c) \ T{c) by contradiction. Assume that there exist 
a plateau of size to > 2 in R[c) \ T{c) i.e. there exists i € i?(c) \ r(c) such 
that Ci — Ci+i = . . . = Ci+m-i- By the hypothesis we know that Q-i > q and 
Ci+m-i > Cm (we assume Cfc+i =0). 

Consider a configuration d such that c G f (d) and c = VJ'{d) for j e [i — 
l,i + m — 1]. Then, d is not LR-decomposable and, by Lemma it does not 
belong to t?(„). A configuration d such that c = Vj{d) for j e [i, i + to] is not 
LR-decomposable either, d ^ Q(^n) ■ The proof for L{c) \ T{c) is very similar. □ 

Lemma 9. Consider n e N and c e G{n)- Let T{c) be the top of c. Any 
LR-decomposition of c is such that both L{c) \ T{c) and R{c) \ T{c) are crazed. 

Proof. Let c be a configuration where R{c) \ T{c) is not crazed, i.e. c contains 
two plateaus not separated by a cliff. Then there are two indices i,j € [t^k], 
i < j, such that Ci_i > q = Q+i, Cj = Cj+i > Cj^2 and for all h £ [i + l,j — 1], 
c/i = Qi+i + 1- Let S = j — i, we prove that c ^ tj(„) by induction on 6. If (5 = 1, 
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Lemma IHI proves the thesis. Suppose the result is true for every S G [l^m], and 
that we have S = m + 1. Let c be a configuration which contains two plateaus 
not separated by a cliff, at distance j ~ i = m + 1. Consider a configuration d 
such that c £ f (d) and c = Vf^{d) for /i e [i — 1, j + 1]. We have the following 
cases: 

• h — i — 1] then di-i > di < di+i, d is not LR-decomposable, d ^ G{n) by 
Lemma 

• h £ — 2]; then there is a plateau in d at position h + 1 and the plateau 
at position j is unchanged, by induction over 6 it holds that d ^ f/(„') ; 

• h = j — I; dj^i > dj < dj+i, d is not LR-decomposable, Lemma [3 says 
that d ^ G{n) ; 

• h — j; there are two plateaus in d at position j — 1 and i, we have that 
d ^ G{n) by induction over 5; 

• h = j + 1; dj^i > dj < dj^i, d is not LR-decomposable, from Lemma[2l 

Choose d such that c G f (d) and c = V^(d) for h G + 2]. The only possible 
values for h are the non-decreasing parts in c, i.e. h £ {i + l,j + 1}: 

• if /i = « + 1; there are two plateaus in d at position i + 1 and j, by induction 
over 6 it holds that d ^ G{n) 

• if /i = j + 1; dj^i > dj < dj+i and from Lemmad d ^ G(n)- 

Therefore among all the ancestors of c which create the plateaus, none is in G{n) , 
hence c ^ G{n) ■ A similar proof can be done if L(c) \ T{c) is not crazed. □ 

It is obvious that the cardinality of T(c) is bigger of equal to 1 for all con- 
figurations. Using very simple examples one can verify that |T(c)| can also be 
equal to 2, 3 or 4. The following result proves that these are the only possible 
values for the cardinality of T{c) when c belongs to an orbit graph. 

Lemma 10. Consider n e N and c e Then |T(c)| < 4. 

Proof. If c = (n) then the thesis trivially holds. Assume that c 7^ (n), then c 
should have an ancestor in G{n)- By contradiction, let \T{c)\ > 4. 

Consider a LR-decomposable configuration d such that c S Vl{d) for i S 
T(c). Then, according to Lemma [3 d has only two possible pre-images 

1) . . . ,d; - + l,d/+i,(i;+2,rf;+3, • • • with = (i;+2 = rf;+3 ; 

2) ■ • ■ , dm-4, dm-3, dm-2, ^m-l + 1, dm — 1, . . . with dm-4 ~ = dm-2 ] 

where I = mmT{d) and m = maxT{d). Both 1) and 2) contradict Lemma|Hl □ 

The following proposition gives a precise characterization of the configura- 
tions of the orbit graph. Its proof is very technical as many different cases have 
to be considered, but each of them is solved quite simply using the previous 
lemmas. 
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Proposition 11. Consider ti e N and c e Q(n)- Then c has a crazed LR-de- 
composition. 

Proof. Let n G N, c = (ci, . . . , Cfc) G Q(„) and T(c) be the top of c. There are 4 
cases, depending on the cardinahty of T(c) fLemma, [T(Hl . 

• If \T{c)\ = 1, then by Lemma|2lany LR-decomposition of c into L(c) and 
i?(c) is crazed (there are no additional plateaus in T(c)). 

• If |T(c)| = 2, choose L(c) = and i?(c) = [t+l,fc] where T(c) = [t,t+l]. 
Since c is LR-decomposable (Lemmad, this is a vaHd LR-decomposition. 
Again there are no plateaus in L{c) T{c) and in R{c) nT{c), hence this 
is a crazed LR-decomposition (Lemma 

• If \T{c)\ = 4, let t e [l,fc] such that T(c) = {t,t + l,t + 2,t + 3}. Choose 
L(c) = + 1] and R{c) = [t + 2,fc]. Again, it is clearly a valid LR- 
decomposition. We prove that both L{c) and R{c) are crazed. The plateau 
T(c) can be obtained from the configurations d G G{n) such that c = Vl'{d), 
i e [t,t4-3]: 

— if i = i or i = t + 1; then di > di+i < (ii+2, which is not possible 
because of Lemma Q 

— ifi = t + 3;dt= dt+i = dt+2 < dt+3, which is not possible because 
of Lemma IHl 

— if i ~ t + 2; dt ~ dt+i < dt+2 and because of Lemma L{d) = 
[1, t + 2] is crazed since L{d) n T{d) — 0. Therefore L{c) is crazed as 
L(c) = L((i) and for all i G i(c), Ci = di. For i?(c), there are two 
cases. 

* Either ct+s > Ct+A+2; then the plateau ct+2,Ct+3 in i?(c)nT(c) is 
separated from any other plateau in R{c) by the cliff at position 
t + 3, hence R{c) is crazed. 

* Or ct+3 = Ct+4 + 1; then dt+2 > dt+3 = dt+A and [t+'3, k] is crazed 
in d. This means that if j > t -I- 4 is the lowest index such that 
dj = dj+i, there is a cliff somewhere at index h, t + A<h<j 
(Lemma El . Hence this cliff is also in c, and there is no plateau 
in c between t + 3 and h. Therefore the plateau ct+2,Ct+3 is 
separated from any other plateau in R{c) by the cliff at index h, 
R[c) is crazed. 

Similar results hold if c = Vl{d), i E [t,t + 3]. Therefore there are only 
two possibilities for d, and for both of them L{c) and i?(c) are crazed. 

• If \T{c)\ = 3, let T{c) ^ {t,t+l,t + 2}. Suppose that X{c) = [t + l,k] is 
crazed, then let L(c) — [l,t] and R{c) = X{c). This is a LR-decomposition 
of c (Lemmad, it is crazed because L(c)nT(c) is also crazed (LemmaEJ. 

If X{c) is not crazed, clearly [t + 2, k] is crazed because of LemmaEl We 
need to prove that Y{c) = [1, t+1] is necessarily crazed. Let j G [t+2, k] be 
the lowest index such that Cj = cj+i . Remark that for all ft- G [t + 2, j — 1] , 
Cfi = Ch+i + 1. If j = t + 2, we are in the case \T{c)\ = 4, solved previously. 
Otherwise, for any ancestor c? of c such that c = V'^(c?), /i G [t,j + 1], it 
holds that 
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— \{ h = t] dt > dt+i < dt+2 which is impossible because of Lemma[3 

— if ft, e [t + 1, j — 1]; dh > d/i+i = (i/i+2, from Lemma El this is not 
possible; 

— if /i = j, dj-i > dj < c?j+i, Lemma 13 proves that it is impossible; 

— ii h = j + 1, plateau at dj-i, dj, by induction (see proof of LemmaEJ 
it leads to the case |T(e)| — 4 for an ancestor e of c, which implies 
that Y{c) is crazed. 

If d is such that c = Vl{d), h € {t,t + l,t + 2, j + 1}, it holds that 

— \{ h = t; dt > dt+i = dt+2, by Lemma El this is impossible; 

— \ih = t+l; the proof is exactly the same as for the sub-case i ^ t + 2 
of the case |T(C)| = 4: Y{c) is crazed; 

-if/i = t-|-2orft = j + l; dh-2 > dh-i < dh which is impossible 
(Lemma CJ. 

Therefore Y{c), X{c) is a valid crazed LR-decomposition of c. 

□ 

The converse of Proposition El is proved using another technical lemma. 

Lemma 12. Consider a configuration c, c ^ (n) for all n G N, which admits 
a crazed LR-decomposition. Then, there exists d such that c G f{{d}) and d 
admits a crazed LR-decomposition. 

Proof. Assume that c = {ci,...,Ck) is such that c ^ (n) for n = J2'i=i'^i- 
Moreover, assume that c admits a crazed LR-decomposition and denote it by 
L(c) and i?(c). 

If L{c) = 0, then if |i?(c)| = 1, nothing can be done: this is the case c = (n), 
which is not possible by hypothesis. Otherwise, build a configuration d as 
follows: di = ci + 1, c?2 = C2 — 1 and di = Ci for i e [3, /c]. Then, c = V{{d) 
and L{d) = 0, R{d) = R{c) is a crazed LR-decomposition of d. Note that if 
C2 = 1, (i2 = so c? has length k—1. In that particular case, we define L{d) ~ 0, 
R{d) = 1]. 

If \L{c) = 1|, we have two cases. If ci > C2, we could have chosen L{c) ^ 9 
and R{c) = case solved previously. Otherwise, define d such that di ~ 

ci — 1, ^2 = C2 + 1 and di = Ci for i e [3, fc]. It holds that c = V2{d), and 
L{d) = L{c), R{d) = R{c) is a crazed LR-decomposition of d. Remark that if 
ci = 1, c?i = so d has to be shifted by 1 to the left. In that particular case 
also, we define L{d) = 0, R{d) = [l,k — 1] which is a crazed LR-decomposition 
of d. 

If \L{c)\ > 1, then we have two cases. 

• L(c) contains at least one plateau: let i be the least index such that 
Ci ~ Ci+i. Define a configuration d as follows: di = Ci — 1, di+i = ci^i + 1 
and dj — Cj for j E [1, k] -\- 1}. Clearly, c = Vl{d) and L{d) ~ L{c), 

R{d) = R[c) is a crazed LR-decomposition by hypothesis. Again, \{ i ~ 1 
and ci = 1, d has to be shifted by 1. Let m = maxL(c), then L{d) = 
L{c) \ {to}, R{d) = {R{c) U {to}) \ {k} is a crazed LR-decomposition of d. 
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• L(c) contains no plateau: define a configuration d such that o?i = ci — 1, 
^2 = C2 + 1 and di = Ci for i G [3, k]. Clearly, c = Vjid), and = L{c), 
R{d) ~ R{c) is a crazed LR-decomposition of d. Once more, if ci = 1, let 
TO = maxL(c), = L(c) \ {to} and R{d) = (i?(c) U {m}) \ {k} for the 
same result. 

□ 

The next proposition proves that having a crazed LR-decomposition is suf- 
ficient to belong to an orbit graph. 

Proposition 13. // a configuration c admits a crazed LR-decomposition, then 
there is a n G N such that c G G{n) ■ 

Proof. If c = (n) for some n G N then we are done. Now, assume that c = 
(ci, . . . , Cfc) is such that c ^ [n) with n = X)i=i Ci- Using LemmaEl build a 
sequence of configurations c?°, d^, . . . , c?'', . . . such that = c, d'* G / for 
/i > and d'' admits a crazed LR-decomposition. Remark that this sequence 
must be finite. Indeed, for all /i G N, by Lemma 01 E{d^'^^) > E{d'^) and, by 
Lemma El E{{n)) > E{d'') if d'' ^ (n). Therefore there is ^ G N such that 
= (n), hence there is a path in G{n) from = (n) to d^ = c. □ 

Because of Proposition any fixed point 11 of G(n) has very precise char- 
acteristics. It admits a crazed LR-decomposition L{U), R{I1), and it has no 
cliffs. Therefore, both L{I1) and i?(n) have at most 1 plateau since they are 
crazed. Moreover, there may be another plateau at the junction between L{Il) 
and i?(n), i.e. at most 3 plateaus in 11. 



The structure of the fixed points is described on Figures O Figure 5(a) 



represent the fixed points 11 such that |r(n)| = 1, Figure 5(b) is for the fixed 
points n such that |r(n)| > 2. 




Figure 5: Structure of the fixed points. 



3.3 A kind of magic 

From Corollary El we know that for any ti G N, the configuration (n) leads to at 
least one fixed point. In this section we compute precisely the number of fixed 
points of SSPM with initial condition (n). 
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In order to understand how a fixed point can be obtained, we try to give a 
visual construction. Consider Figures |3 The n grains of the fixed point must 
be arranged in the grayed part and can partially occupy the dashed frame with 
the supplementary constraint that grains in the dashed part must be as much 
clus tered to the ground as possible. Boxes labeled 1, 2 and 3 in Figures 5(a) 



and 5(b) cannot be filled (for more details see the proof of Lemma^j). Remark 
that if p is the height of the grayed part, then this area contains grains in 
Figure [5'(a)| and +p grains in Figure 5(b) Appendix El shows all the possible 
fixed points reachable from the initial condition (n) for n £ [1,32]. 

Lemma El will be the main tool that we use to count the number of fixed 



points, it proves that the "shapes" outlined in Figure 5(b) describe exactly all 
the possible fixed points. Let gi{n) be the numbers of fixed points 11 such that 
|r(n)| = 1, and g2(n) the numbers of fixed points 11 such that |r(n)| > 2. 

Lemma 14. For any n G N, consider SSPM with initial condition {n). The 
number of fixed points of is given by G{n) ~ gi{n) + g2{n). 

Proof. Let n g N, and 11 be a fixed point of (n). If |r(n)| = 1, clearly 11 can be 
constructed as shown in Figure 5(a) (at most one plateau on the left and one 
on the right). Moreover, it cannot be constructed from Figure 5(b) since the 

boxes at the top (labelled 1 and 2) are left empty. 

If |T(n)| > 2, the fixed point is not represented on Figure 5(a) since the 
boxes labelled 1 and 2 cannot be filled. To show that it is constructible from 
Figure 5(b)| let L(I{) = [l,t] and i?(n) = [t + 1, fc] be a crazed decomposition 
of n = (Hi, . . . , Ilfc). If lit = n(_|_i, cut Figure 5(b) in two parts at the middle 
of the configuration. It is clear that Liji) fits in the left part (at most one 
plateau), and i?(n) fits in the ri ght part for the same reason. If XIj = IIj+i + 1, 
T(n) C i(n). Cut Figure [5(b)| in two parts, at the right of the two grains on 



top of the grayed pile, L{J\) and Riji) fit. The symmetrical case is similar. 

Conversely, all the configurations with n grains constructible from Figures 
are clearly fixed points, are LR-decomposable, and hence are fixed points of (n) 
(Proposition E|l ■ Therefore, the total n umbe r of fixed points is the number of 
configurations constructib le fro m Figure 5(a) plus the number of configurations 
constructible from Figure 5(b) i.e. G(n) ~ gi{n) + g2{n). 



□ 



The two following lemmas give the exact expression of gi{n) and g2{n). 

Lemma 15. For any n e N, consider SSPM with initial condition (n). The 
number of fixed points of Q(^n) with top of length 1 is given by 

^ 1 ifn — p^<p—l, 



p- 



9i{n) 



p^ — 1 ifp <n — p^ < 2p — 1 
otherwise, 



where p is the unique integer such that p^ < n < {p + 1)^. 

Proof. For n g N, consider a fixed point 11 g G{n)- Since by hypothesis |T(n) 



1|, the overall structure is illustrated in Figure 5(a) As 



fixed, we can 

determine p: it is the unique integer satisfying p^ < n < {p+ 1)^, i.e. p = [\/?iJ • 
Now, let u = n—p^ be the number of grains left after having arranged the grayed 
zone. Distributing these u grains consecutively and in all possible manners on 
the borders of the grayed zone starting from bottom to top gives all possible 
fixed points (lemma ^J. To be more precise we must distinguish three cases: 
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• < u < p: we put all u grains in the free boxes on the left, from bottom 
to top; this gives a fixed point. Then we put only u — 1 grains on the left 
and 1 in the free box at the bottom on the right; this gives another fixed 
point. This process is iterated until there are gains on the left and u on 
the right. It is clear that this procedure gives u + 1 fixed points. 

• p < u < 2p — 2: we start by putting p — 1 grains in the free boxes on the 
left, and the remaining u — p + 1 grains in the free boxes on the right, 
starting from bottom to top; this gives a fixed point. Then we put only 
p—2 grains on the left and the u — p + 2 remaining grains in the free boxes 
on the right, proceeding from bottom to top; this gives another fixed point. 
Then, we can put p — 3 grains on the left and so on until there are p — 1 
grains on the right. It is clear that {p ~ 1) — {u — p + 1) + 1 = 2p — u — 1 
distinct fixed points can be generated in this manner. 

• u = 2p — 1 or u = 2p; then there are necessarily p grains on the left or on 



the right, hence there is one grain in box 1 or 2 (Figure 5(a) I. This should 
not be allowed, as it would mean that |T(n)| > 1, which is not the case. 
In this last case, there are fixed points. 

Finally, remark that box number 3 is not taken into account either, since it 
would mean that all dashed boxes are filled and therefore we would have chosen 
p + 1 as height of the pile instead of p. □ 

Lemma 16. For any n e N, consider SSPM with initial condition (n). The 
number of fixed points of Q(n) with top of length bigger than 1 is given by 

■ — p^ — p + I ifn—p^—p<p—l , 
92{n)^{p ifn-p^-p = p, 

n + p^ + 1 ifp +l<n — p^— p<2p+l , 

where p is the unique integer such that p^ + p <n < (p + 1)^ + (p + 1). 

Proof of \lfi[ The proof is similar to the one of Lemma^J For n e N, consider a 
fixed point c g G{n)- By hypothesis \T{c)\ > 2, therefore the overall structure is 




the one illustrated in Figure 5(b) Since n is fixed, we can determine p: it is the 
unique integer satisfying p^ +p < n < (p+ 1)^ + (p+ 1). Let -u = n—p^ —p be the 
number of grains left after having arranged the grayed zone. Distributing the v 
grains in all possible ways gives the number of fixed points 11 with |r(n)| > 2 
(LemmaEJ- Again, we must distinguish three cases: 

• < w < p: we put all v grains in the free boxes on the left from bottom 
to top; this gives a fixed point. Then we put only v — 1 grains on the left 
and 1 in the free box on the right at the bottom; this gives another fixed 
point. This can be iterated until there are grains on the left and v on 
the right. It is clear that this procedure gives v + 1 fixed points. 

• V = p: again, we can put p grains on the left, on the right and so on 
until there are on the left, p on the right. Therefore there should be 
p + 1 fixed points, but in fact the first one and the last one are exactly 
the same: top of length 3, and no plateaus. What happens is that the 
reference column is not at the same position, it is shifted by one, which 
does not matter. This is the only case of duplicated fixed point. 
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• p < V < 2p + 1: we start by putting p grains in the free boxes on the left, 
and the remaining v — p grains in the free boxes on the right, starting from 
bottom to top; this gives a fixed point. Then we put only p — 1 grains on 
the left and the v — p+1 remaining grains on the free boxes to the right, 
proceeding from bottom to top; this gives another fixed point. We proceed 
until there arep grains on the right, it is clear that p—{v—p) + l = 2p—v+l 
distinct fixed points can be generated in this manner. 

Remark that boxes at the top of the dashed pile must not be taken into account 
in the computation of the number of fixed points for it would mean that all 
dashed boxes are filled and therefore we would have chosen p + 1 as height of 
the pile and not p. Moreover, if only one of them is filled, we generate a fixed 
point already counted in Lemma, ITfil (ton of length 1). □ 

The following proposition gives a closed formula for the number of fixed 
points in the orbit of initial condition (n) . The formula is somewhat "magical" 
since it is very simple but we have neither practical nor visual explanation for 
it. 

Proposition 17. For any n G N, consider SSPM with initial condition (n). 
The number of fixed points of G{n) is given by G{n) ~ Lv^J ■ 

Proof. Let n S IM and p = [VnJ the only integer such that p^ < n < {p + 1)"^ , 
there are three cases which correspond to the cases of Lemmas El and El 

• li p^ < n < p'^ + p — 1; then (p — 1)^ + (p — 1) < n < p'^ + p. This is 
case 1 for Lemma El and case 3 for Lemma El hence from Lemma El 

G{n) = [n -p'^ + l] + [3(39 - 1) - n + {p - if + I] = p = [^^J . 

• If p^ + p < n < p'^ + 2p — I; we are in case 2 of Lemma El and case 1 of 
LemmaEl hence from Lemma ll4[ G{n) = [2p— n+p^ — + p^— p+1] = 
P = [Vn\ . 

• If n = p^ + 2p, this is case 3 for Lemma[T5land case 2 for LemmaEl From 
Lemma ITU we find G{n) — 0+p — [^/n\ . 

□ 

Remark that it would also be possible to give the exact expression of each 
of these fixed points, but it would be complex and of no interest here. To have 
an idea of what they look Hke, please refer to Appendix El 

Finally, remark that we did not take into account the initial position of the 
columns. For the same fixed point, there may exist different fixed points which 
have the same shape, but at different indices. In this paper we do not consider 
this fact, we only take into account the general shape of the configurations. 

4 Conclusions and future work 

In this paper we have introduced SSPM: a symmetric version of the well-known 
SPM model. We have proved that SSPM has fixed point dynamics and exhibited 
the precise structure of the fixed points which are in the orbit of initial condition 
(n). Moreover, we showed a simple closed formula for counting the number of 
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distinct {i.e. having different shape) fixed points. Remark that this result is 
surprising since the combinatorial complexity of the orbit graphs becomes higher 
and higher when the number n of grains grows. This complexity contrasts with 
the simplicity of the formula for the number of fixed points: l\/n\ . Moreover, 
this formula is to some extent fascinating: although it is very simple, we have 
neither a practical nor a visual explanation for it. 

This research can be continued along three main directions: 

• Corollary says that, starting from any initial configuration, SSPM has 
fixed point dynamics. Can we give a formula or at least tight bounds for 
the shortest path to a fixed point? For the longest? 

• Section 13.21 gives a precise characterization of orbit graphs for initial con- 
ditions made of one single column. It would be interesting to extend this 
characterization to more general initial conditions or at least to find an 
alternative characterization. 

• The model we introduced is intrinsically sequential: only one grain moves 
at each time step. It would be interesting to introduce a model similar to 
SSPM but with synchronous update. This would be even more realistic 
than SSPM for the simulation of natural phenomena. 
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Fixed points of (n) for 1 < n < 32 

n = 1 n = 2 n = 3 n = 4 
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n = 29 n = 30 n = 31 n = 32 
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